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ABSTRACT 

Aims. The chemical evolution of galaxies on a cosmological timescale is still a matter of debate despite the increasing number of 
available data provided by spectroscopic surveys of star-forming galaxies at different redshifts. The fundamental relations involving 
metalUcity, such as the mass-metalUcity relation (MZR) or the fundamental-metallicity relation, give controversial results about the 
reality of evolution of the chemical content of galaxies at a given stellar mass. In this work we shed some light on this issue using 
the completeness reached by the 20k bright sample of the zCOSMOS survey and using for the first time the nitrogen-to-oxygen ratio 
(N/O) as a tracer of the gas phase chemical evolution of galaxies that is independent of the star-formation rate. 
Methods. Emission-line galaxies both in the SDSS and 20k zCOSMOS bright survey were used to study the evolution from the local 
Universe of the MZR up to a redshift of ~ 1.32, and the relation between stellar mass and N/O (MNOR) up to a redshift of ~ 0.42 
using the N2S2 parameter. All the physical properties derived from stellar continuum and gas emission-lines, including stellar mass, 
star formation rates, metallicity and N/O, were calculated in a self-consistent way over the full redshift range. 
Results. We confirm the trend to find lower metaUicities in galaxies of a given stellar mass in a younger Universe. This trend is even 
observed when taking possible efl'ects into account that are due to the observed larger median star formation rates for galaxies at 
higher redshifts. We also find a significant evolution of the MNOR up to z ~ 0.4. Taking the slope of the O/H vs. N/O relation into 
account for the secondary-nitrogen production regime, the observed evolution of the MNOR is consistent with the trends found for 
both the MZR and its equivalent relation using new expressions to reduce its dependence on star-formation rate. 
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1. Introduction 

Chemical abundances provide important clues to the evolution- 
ary history of galaxies along cosmic time since they are the re- 
sult of the joint action of several physical processes such as su- 
pernova feedback, gas inflow/outflow, mergers, and interactions, 
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Even if it is not easy to disentagle the effect of each individual 
process, it may be possible to constrain galaxy formation theo- 
ries by studying how chemical abundances for galaxies of dif- 
ferent masses and/or in different environments evolve as a func- 
tion of cosmic epoch. Indeed, as cosmological time progresses, 
galaxy evolution models predict that both the mean metallic- 
ity and stellar mass of galaxies increase with age as galaxies 
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undergo chemical enrichment and grow through merging pro- 
cesses. At any given epoch, the accumulated history of star for- 
mation, gas inflows, and outflows, affects a galaxy mass and its 
metallicity. One therefore expects these quantities to be corre- 
lated in some way and this correlation to provide crucial infor- 
mation about the physical processes that govern galaxy forma- 
tion. 

First discovered for dwarf irregular galaxies (Lequeux et al. 
1979), the relation between the mass (or luminosity) and the 
metallicity for different galaxy populations has been the subject 
of numerous studies (e.g. Skillman et al. 1989; Brodie & Huchra 
1991;Zaritskyetal. 1994; Garnettetal. 1997; Pilyugin & Ferrini 
2000; Contini et al. 2002). The mass-metallicity relation (here- 
after MZR) is now well established in the local universe, thanks 
to the works of Tremonti et al. (2004) based on SDSS and 
Lamareille et al. (2004) based on 2dFGRS data, showing that 
the metallicity of galaxies tends to increase with their mass or 
luminosity. This trend has been shown to extend to much lower 
galaxy masses (Lee et al. 2006; Saviane et al. 2008), confirming 
the idea that a single mechanism may govern galaxy metallicities 
across five order of magnitude in stellar mass. However, a spe- 
cific environment and/or high star formation rate can explain the 
observed off'set of galaxy population with respect to the global 
MZR. For instance, nuclear inflows of metal-poor interstellar gas 
triggered by galaxy interactions can account for the systemati- 
cally lower central oxygen abundances observed in interacting 
galaxies (Perez et al. 2011; Torrey et al. 2012). 

Hierarchical galaxy formation models that take the chemical 
evolution and feedback processes into account are able to re- 
produce the observed MZR in the local universe (e.g. De Lucia 
et al. 2004; De Rossi et al. 2007; Finlator et al. 2007; Dave 
et al. 2012). However these models rely on free parameters, 
such as feedback efficiency, which are not yet well constrained 
by observations. Alternative scenarios proposed to explain the 
MZR include low star formation efficiency in low-mass galaxies 
caused by supernova feedback (Brooks et al. 2007) and a vari- 
able stellar initial mass function that is more top-heavy in galax- 
ies with higher star formation rates, thereby producing higher 
metal yields (Koppen et al. 2007). 

The evolution of the MZR on cosmological timescales is 
now predicted by semi-analytic models of galaxy formation, 
which include chemical hydrodynamic simulations within the 
standard A-CDM framework (De Lucia et al. 2004; Dave & 
Oppenheimer 2007; Sakstein et al. 201 1). Reliable observational 
estimates of the MZR of galaxies at different epochs (hence 
different redshifts) may thus provide important constraints on 
galaxy evolution scenarios. Estimates of the mass-metalUcity - 
or luminosity-metallicity - relation of galaxies up to z ~ 1.5 
have been already derived but have been limited, until recently, 
to small samples (e.g. Kobulnicky et al. 2003; Liang et al. 2004, 
Maier et al. 2004, 2005, 2006; Hammer et al. 2005; Savaglio et 
al. 2005; Lamareille et al. 2006; Liu et al. 2008; Queyrel et al. 
2009; Yabe et al. 2012). Recent studies of the MZR (Lamareille 
et al. 2009; Perez-Montero et al. 2009; Cowie & Barger 2008; 
Moustakas et al. 2011; Zahid et al. 2011) have been performed 
on larger samples (> 1000 galaxies) thanks to the large and 
deep spectroscopic surveys (VVDS, DEEP2, GOODS, AGES, 
etc) devoted to the formation and evolution of galaxies on cos- 
mological timescales. The general observational result of these 
studies is that the MZR evolves with redshift, in the sense that, 
on average and for a given stellar mass, high-redshift galaxies 
are characterised by lower metallicities. Whether the MZR also 
evolves in terms of shape/slope is still matter of debate. Indeed, 
SavagUo et al. (2005) conclude that there is a steeper slope in 



the distant universe, whereas Lamareille et al. (2009) and Perez- 
Montero et al. (2009) found flattening of tiie MZR at redshifts 
0.7 < z < 1.0. At higher redshifts Erb et al. (2006) derived a 
MZR at z ~ 2 lowered by 0.3 dex in metallicity compared with 
the local estimate, a trend that could extend up to z ~ 3-4 
(MaioUno et al. 2008; Mannucci et al. 2009). 

However, the redshift evolution of the gas-phase metalUc- 
ity in galaxies has been questioned recently by Mannucci et 
al. (2010) and Lara-Lopez et al. (2010). They discovered that 
metalhcity depends not only on the stellar mass, but also on the 
star formation rate (SFR): for a given stellar mass, galaxies with 
higher SFR systematically show lower metallicities. This is the 
so-called fundamental metallicity relation (FMR), i.e., a tight 
relation between stellar mass, gas-phase metallicity, and SFR. 
Local SDSS galaxies show very small residuals (~ 0.05 dex) 
around this relation (Yates et al. 2012; Brisbin & Harwit 2012). 
According to Mannucci et al. (2010), the FMR does not evolve 
with redshift up to z ~ 2.5. This result was later confirmed by 
Cresci et al. (2012) by studying the zCOSMOS 10k bright sam- 
ple up to a redshift » 0.8. This would suggest that the observed 
evolution of the MZR is due to selection effects and to the well- 
established increase in the average SFR with redshift (Noeske et 
al. 2007; Elbaz et al. 2007; Daddi et al. 2007). 

The study of other indicators of the chemical content of 
galaxies and their relation with steUar mass and SFR have not 
been sufficiently explored so far so can shed some light on the 
issue of the evolution of metallicity with cosmic age. This is the 
case of the nitrogen-to-oxygen abundance ratio (N/O), which 
offers some clear advantages over studies of the oxygen abun- 
dance alone. On one hand, N is produced mainly by low- and 
intermediate-mass stars (Henry et al. 2000) in contrast to O, 
which is produced only by massive stars. This makes N/O a 
suitable indicator for the chemical evolution of single starbursts 
(Edmunds & Pagel, 1978; Pilyugin et al. 2003). On the other 
hand, since N is a secondary element {i.e. its yield depends on 
the previous amount of carbon and oxygen in the stars) for the 
majority of the metallicity range, its relation with a primary ele- 
ment, such as O, is relatively independent of the chemodynami- 
cal effects, such as outflows or inflows (Edmunds, 1990, Koppen 
& Hensler, 2005). These effects are the basis of the relation be- 
tween metallicity and SFR, so the relation between N/O and stel- 
lar mass (hereafter MNOR) can be fundamental for studying the 
chemical evolution of galaxies without worrying about the se- 
lection effects of the sampled galaxies at high z. The MNOR for 
galaxies in the local universe has already been investigated by 
Perez-Montero & Contini (2009, hereafter PMC09) who found 
an increase in N/O with stellar mass, since most of the sampled 
SDSS galaxies lie in the metalUcity range for the production of 
secondary N. This relation was used by Amorin et al. (2010) 
to study the chemical evolution of compact galaxies with very 
high specific star formation rates (the green pea galaxies). Later, 
Thuan et al. (2010) explored the evolution of the relation be- 
tween stellar mass and N/H up to z = 0.3, exploiting the fact that 
the production rate of N is faster than for O and its evolution is 
thus less sensitive to calibration uncertainties. 

In this work we investigate both the MZR with its relation to 
SFR and the MNOR in the Local Universe and their evolutions 
with cosmic age. This study is based on the SDSS sample for the 
local universe and the zCOSMOS 20k bright sample for high 
redshifts (up to a z ~ 1.32 for the MZR and z « 0.42 for the 
MNOR). 

This paper is organised as follows. In section 2 we describe 
the SDSS and zCOSMOS samples of the star-forming galax- 
ies used to perform our study of the evolution of both the MZR 
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and the MNOR. The selection of star-forming galaxies and the 
derivation of physical properties, including stellar mass, star for- 
mation rate, oxygen abundance, and the N/O abundance ratio 
are described in section 3. In section 4 we give our results and 
discuss the evolution with cosmic age of the MZR, with a cor- 
rection to take selection effects due to SFR into account, and of 
the MNOR. Finally, in section 5 we summarise our work and 
give our conclusions. 

Throughout this paper we normalise the derived stellar 
masses and the absolute magnitudes with the standard A-CDM 
cosmology, i.e., h - 0.7, Q.,„ - 0.3 and Oa = 0.7 (Spergel et al., 
2003). 

2. Sample selection 

2. 1 . The parent SDSS local sample 

Ss reference sample for the local universe, we used the Data 
Release 7 of the Sloan Digital Sky Survey (hereafter SDSS), 
taking the emission-line measurements from the Max Planck 
Institute for Astrophysics-Johns Hopkins University (MPA- 
JHU) catalogutj] The emission-line measurement of this sam- 
ple is described in Brinchmann et al. (2004). From the whole 
sample of galaxies, we removed duplicated objects and all those 
emission-line galaxies whose H/3, [On] AA 3727, 3729, [Oiii] 
AA 4959, 5007, Ha, [Nii] A 6584, and [S ii] AA 6717,6731 have 
signal-to-noise ratio lower than 2. The choice for this low S/N 
was motivated by the search for a coincidence between the mass- 
metallicity relation of the SDSS sample derived in this work and 
the zCOSMOS subsample at low z, avoiding biases in the sam- 
ple selection for the derivation of metallicity, as discussed be- 
low. The selected sample comprises a total of 299 479 galaxies 
in the 0.02 < z < 0.20 range, but with a mean value z = 0.068. 
The lower redshift limit is caused by the lower limit of the spec- 
tral coverage of the SDSS sample (3800 A), which prevents the 
detection of the [On] A 3727 emission-line at z < 0.02. This 
lower limit also prevents the sample from being contaminated 
by nearby objects with an inaccurate determination of their inte- 
grated properties due to the limited size of the fibre in the SDSS. 
as they were not totally covered by the 3" fibre of SDSS. All 
collected emission lines were extinction-corrected using the law 
by Cardelli et al. (1989) and the reddening constants from the 
Balmer decrement between Ha and H/3 as compared to the the- 
oretical value for mean nebular conditions given by Storey & 
Hummer (1995). 



2.2. The parent zCOSMOS 20k sample 

The COSMOS survey is a large HST-ACS survey, with /-band 
exposures down to /ab = 28 on a field of ~ 2 deg^ (Scoville et 
al. 2007). The COSMOS field has been the object of extensive 
multiwavelength ground- and space-based observations span- 
ning the entire spectrum: X-ray, UV, optical/NIR, mid-infrared, 
mm/submillimetre and radio, providing fluxes measured over 
30 photometric bands (Hasinger et al. 2007; Taniguchi et al. 
2007; Capak et al. 2007, Lilly et al. 2007; Sanders et al. 2007; 
Bertoldi et al. 2007; Schinnerer et al. 2007; Koekemoer et al. 
2007; McCracken et al. 2010). 

The main spectroscopic follow-up survey undertaken in the 
COSMOS field, zCOSMOS (Lilly et al. 2007), used 600 hours 
of ESO observing time with the VIMOS multi-object spec- 
trograph (Le Fevre et al. 2003) mounted on the Melipal 8m- 
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telescope of the ESO-VLT. The zCOSMOS spectroscopic survey 
consists of two parts: zCOSMOS-bright and zCOSMOS-deep. 
The zCOSMOS-deep targets ~ 10000 galaxies within the cen- 
tral 1 deg- of the COSMOS field, selected through colour cri- 
teria to have 1.4 ^ z ^ 3.0. The zCOSMOS-bright is purely 
magnitude-limited and covers the whole area of 1.7 deg^ of 
the COSMOS field. It provides redshifts for ~ 20000 galaxies 
down to /ab ^ 22.5 as measured from the HST-ACS imaging. 
The success rate in redshift measurements is very high, 95% in 
the redshift range 0.5 < z < 0.8, and the velocity accuracy is 
~ 100 km/s (Lilly et al. 2009). Each observed object has been 
assigned a flag according to the reliability of its measured red- 
shift. Classes 3.x, 4.x redshifts, plus Classes 1.5, 2.4, 2.5, 9.3, 
and 9.5 are considered a secure set, with an overall reUabiUty of 
99% (see for details Lilly et al. 2009). 

For this work, we used the zCOSMOS-bright survey final re- 
lease: the so-called 20k sample for about 20000 galaxies with 
z < 2 and secure redshifts according to the above flag clas- 
sification (18 206 objects in total, irrespective of redshift and 
including stars). Observations for the zCOSMOS-bright survey 
were acquired with the medium-resolution {R = 600) grism of 
VIMOS, providing spectra over the 5550 - 9650A wavelength 
range. The total integration time was set to one hour to secure 
redshifts with a high success rate. The spectral range of observa- 
tions enables following important diagnostic emission lines to 
compute metallicity up to redshift z ~ 1.5. The observations 
were acquired with a seeing lower that 1.2". 

2.2.1. Emission-line measurement 

Spectroscopic measurements (emission and absorption lines 
strength) in zCOSMOS were performed through the automated 
pipeline platefit-vimos (Lamareille et al., in preparation) similar 
to the one performed on SDSS {e.g. Tremonti et al. 2004) and 
VVDS spectra (Lamareille et al. 2009). This routine processes 
spectra in two steps. The stellar component of galaxy spectra is 
fitted as a combination of 30 single stellar population templates, 
with different ages and metallicities, from the library of Bruzual 
& Chariot (2003) resampled to the velocity resolution of zCOS- 
MOS spectra. After removal of the stellar component, emission 
lines are fitted together as a single nebular spectrum made of a 
sum of Gaussians at specified wavelengths. Further details can 
be found in Lamareille et al. (2006, 2009). 

For this work we selected those zCOSMOS emission-line 
galaxies with an S/N higher than two for the involved emission- 
lines in each redshift regime, which leaves a total of 5 331 ob- 
jects. 



3. Physical properties 

3. 1 . Selection of star-forming galaxies 

We selected star-forming galaxies from both the parent SDSS 
and 20k zCOSMOS samples by means of several empirical di- 
agnostic diagrams based on bright emission-line ratios. We iden- 
tified narrow-line AGNs (Seyfert 2 and LINERs), taking into ac- 
count that the broad-line AGNs were already removed. The diag- 
nostic diagrams used for the final classification are not the same 
for all the galaxies, as the available emission lines depend on the 
considered redshift range. Indeed, the strategy of the zCOSMOS 
survey, with a fixed 5550-9650 A wavelength range, implies that 
the set of bright emission-lines used in the diagnostic diagrams 
changes with redshift. This is also true for the oxygen abundance 
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Fig. 1. Diagnostic diagrams used to select star-forming objects in the studied sample in each redshift range. From left to right and 
from up to down: [N ii]/Ha vs. [S n]/Ha, [N ii]/Ha vs. [O ni]/H/3, [O ii]/H/3 vs. [O m]/Hy6, and [O n]fRf3 vs. [Ne in]Hy6. 



determination, which is based on the same sets of emission lines 
(see section [34] l. 

For the 315 zCOSMOS emission-line galaxies in the redshift 
range 0.01 < z < 0.13, we used the Ha classification proposed 
by Lamareille (2007) based on [Nii], [Sii], and Ha emission- 
line ratios and independent of the reddening. According to this 
classification, objects are classified as narrow-line AGNs if 

log([Sii]/Ha) > -0.3, and 

log([Nn]/Ha) > -0.4 (1) 
or 

log([Nii/Ha') > -1.05 ■ log([S ii]/Ha). (2) 

This gives a total of 19 AGNs (6%), which are rejected from the 
sample to be analysed. The [N iiJ/Ho' vs. [S ii]/Hq; plot for the 
sample of zCOSMOS 20k emission-line galaxies in this redshift 
regime is shown in upper left hand panel of Fig. [T] 



For both the redshift range 0.15 < z < 0.44 in the zCOS- 
MOS sample and for the SDSS parent sample, we used the well 
known and commonly used diagnostic diagrams (e.g. Baldwin 
et al. 1981; Veilleux & Osterbrock 1987) based on the line ratios 
[O iii]/Hy6 and [S ii]/Ha. In this diagnostic diagram, all points be- 
low the separation curve defined by Kewley et al. (2001) can 
be considered as purely star-forming objects. Since our analysis 
for both oxygen and nitrogen abundances are based mainly on 
[N n] emission-line we preferred not to use the diagram that de- 
pends on the [N iil/Ha for two reasons: i) this diagnostic imposes 
an upper limit on the N2 = log([Nii]/HQ') parameter, and hence 
for the upper metallicity and ii) according to PMC09, both the 
composite (i.e a region in the diagram where are objects with a 
mixture of ionising sources, Kewley et al. 2006) and the AGN 
region can be contaminated with pure star- forming objects with 
very high N/O. The total number of zCOSMOS galaxies classi- 
fied as narrow-line AGN in this redshift range is 188, which is 
9.6% of the sample. The right 
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hand upper panel of Fig. [T] shows the sample, along with the 
curve of separation. The limit log([Oni]/HjS) = 0.5 shows the 
separation between Sey2 and LINER galaxies. Regarding SDSS, 
a total of 44729 objects were classified as narrow-line AGNs 
(17.5 % of the sample) and 254 750 objects were selected as star- 
forming galaxies for the analysis of the fundamental relations in 
the local universe. 

At higher redshift (0.47 < z < 0.9), the Ha, [Nii] A 6584, 
and [Sii] AA 6717,6731 emission lines are not visible anymore 
in the observed wavelength range of zCOSMOS spectra. In this 
redshift range, we instead used the blue diagnostic diagrams, as 
defined by Lamareille et al. (2004), involving the [Oiii]/HyS and 
[Oii]/HjS emission-line ratios. Unlike the diagnostic diagrams 
used for lower redshift, this method does depend on reddening 
due to the long wavelength baseline between [O ii] A 3727 A and 
Hy6. One way to minimise this dependence is to use the ratio of 
equivalent widths instead of fluxes as proposed by Kobulnicky et 
al. (2003), but this method must be corrected owing to the differ- 
ent shape of the underlying stellar continuum in each galaxy. To 
do so, we took the expression proposed by Liang et al. (2007), 
who propose the use of a multiplicative factor a to the ratio of 
equivalent widths of [On] and H/? as a function of the spec- 
tral stellar break at 4000 A. In our case, this break was auto- 
matically measured in all the zCOSMOS 20k spectra at this z 



Table 1. Relative number of star-forming galaxies and narrow- 
line AGNs in each redshift bin in the zCOSMOS bright sample. 



by platefit-vimos in the same way as described in section 2.2.1 
By considering the separation curve proposed by Marocco et al. 
(2011), plotted in the left hand lower panel of Fig. [T] we clas- 
sified 262 galaxies (9.9%) as narrow-line AGNs, and thus we 
rejected them. 

Finally, for the highest redshift galaxies {z > 0.9), we used 
the method defined by Perez-Montero et al. (2007) and based on 
the [On] A 3727 and [Nem] A 3869 emission lines in relation 
to the brightest available Balmer hydrogen emission line. This 
diagnostic also utilizes the empirical relation between the emis- 
sion line fluxes of [O iii] A 5007 and the [Ne iii] line at a bluer 
wavelength for a sample of well-characterised ionised gaseous 
nebulae of the Local Volume. The ratio found by Perez-Montero 
et al. (2007) [F([Oin])/F([Neiii]) ^ 15.37] agrees within the er- 
rors with the ratio measured in those 20k zCOSMOS galaxies 
with a trustable measurement of the two lines. For those galax- 
ies with no measurement of H/3, we used Hy or H6 and the the- 
oretical relations between Balmer lines from Storey & Hummer 
(1995) coefficients. We also used this method for those objects 
in the range 0.5 < z < 0.9 with no confident measurement of the 
[O III] emission lines. Regarding extinction correction, unlike the 
previous redshift range where the ratio of equivalent widths was 
used, here this does not work due to the stellar continuum vari- 
ations at [Ne III] wavelength, according to Perez-Montero et al. 
(2007). Instead, we used the Balmer decrement in those galaxies 
with more than one trustable hydrogen emission-line and, other- 
wise, we took the derived inner stellar extinction from the spec- 
tral synthesis fitting. We also took the separation curve proposed 
by Marocco et al. (201 1), as shown in the lower right-hand panel 
of Fig.[T] The number of rejected objects in this regime is much 
higher than in the other three (190, or 46.0% of the sample), but 
this cannot be taken as a much higher fraction of AGNs in this 
redshift range. On the contrary, it is probably a selection effect 
as AGNs have higher S/N for the involved lines than pure star- 
forming galaxies. In fact, increasing the S/N cut to 3, 5, and 10 
increases the fraction of AGNs to 53%, 64%, and 80%, respec- 
tively. 



z bin 


Total 


SF 


NL-AGN 


(%) 


0.0-0.2 


485 


460 


25 


5.2 


0.2 - 0.4 


1637 


1487 


150 


9.2 


0.4 - 0.6 


766 


640 


126 


16.5 


0.6-0.8 


1799 


1590 


209 


11.6 


0.8-1.0 


568 


445 


123 


21.7 


1.0-1.32 


76 


46 


30 


39.5 


All 


5331 


4668 


663 


12.4 



3.2. Stellar mass 




0.2 0.4 0.6 0.8 1 

Redshift 



1.2 



Fig. 2. Relation between redshift and stellar mass (in solar 
masses) for the complete zCOSMOS 20K sample (brown points) 
and star-forming selected galaxies (violet points). The solid lines 
show the logarithmic fitting to the limiting masses of the star- 
forming sample for levels 25%, 50%, and 75% of completeness. 
Lower and left hand panels also show the distributions for both 
the complete sample (empty histogram) and the star-forming 
galaxies (filled histogram) of z and M*, respectively. 



Total stellar masses for galaxies in the SDSS parent sam- 
ple were taken from the MPA/JHU catalogue, based on fits to 
the photometric SED. Regarding masses in the zCOSMOS sam- 
ple, they were derived from fitting stellar population synthesis 
models to both the broad-band optical/near-infrared (CFHT: u, 
i, Ks; Subaru: B, V, g, r, i, z; Capak et al. 2007) and far-infrared 
(Spitzer/IRAC: 3.6yum, 4.5jum; Sanders et al. 2007) photome- 
try, and two spectral features (H5 absorption line and D(4000) 
break) when observed in the VIMOS spectra, using a chi-square 
minimisation for each galaxy. The different methods used to 
compute stellar masses, based on diff'erent assumptions about 
the population synthesis models and the star formation histo- 
ries, are described in Bolzonella et al. (2010). The accuracy of 
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Table 2. Derived minimum masses (in log(M/MQ)) for different 
redshift bins obtained for different levels of completeness 



z bin 


M;„,„ (25%) 


ML„ (50%) 


M;„„ (75%) 


0.0 - 0.2 


7.45 


7.63 


7.80 


0.2 - 0.4 


8.49 


8.70 


8.96 


0.4 - 0.6 


8.99 


9.17 


9.49 


0.6-0.8 


9.28 


9.51 


9.85 


0.8-1.0 


9.52 


9.76 


10.11 


1.0-1.32 


9.71 


9.95 


10.32 



the photometric stellar masses is satisfactory overall, with typ- 
ical dispersions due to statistical uncertainties and degenera- 
cies of the order of 0.2 dex. The addition of secondary bursts 
to a continuous star formation history produces systematically 
higher (up to 40%) stellar masses, while population synthesis 
models taking the TP-AGB stellar phase into account (Maraston 
2005) systematicaly lower M* by 0.10 dex. Finally, the uncer- 
tainty on the absolute value of M* due to assumptions on the 
Initial Mass Function (IMF) is within a factor of 2 for the typical 
IMFs usually adopted. In this paper, we have adopted the stellar 
masses calculated with the stellar population models of Bruzual 
& Chariot (2003), with the addition of secondary bursts to the 
standard declining exponantial star formation history. 

The left hand panel of Fig. 2 shows the stellar mass distribu- 
tion (in units of solar mass) of the complete zCOSMOS 20k sam- 
ple of 18 206 galaxies with secure redshifts. The median value of 
this distribution is log(M*) - 9.96. The stellar mass distribution 
for the 4 668 star-forming selected galaxies has a median value 
of 9.66. The lower panel in the same figure shows the redshift 
distributions for the same samples, with mean values of 0.59 
for the complete sample and 0.52 for the selected star-forming 
galaxies. 

For the analysis of stellar masses, we must consider that the 
minimum "detected" stellar mass is a function of redshift. We 
thus defined a minimum mass for each redshift bin, which is the 
lowest mass at which the mass function can be considered reli- 
able and unaffected by incompleteness on mass-to-light (Ilbert et 
al. 2004, Pozzetti et al. 2007). To calculate the minimum mass, 
we first computed the limiting mass, which is the stellar mass 
that an object would have at the limiting magnitude, following 
the expression 

log(ML) = log(M*) + 0.4 ■ (Us - he,) (3) 

where M,*^^ is the limiting mass in solar masses, M* the stellar 
mass in solar masses, lobs the observed /-band magnitude, and 
Isei the /-band magnitude limit for each field, which is equal to 
22.5 for the studied zCOSMOS sample. In Fig.|2]we show the 
stellar masses of the selected star-forming galaxies as a function 
of redshift in the range < z < 1 .4. In the same plot we show 
the logarithmic fits to the 25, 50, and 75 percentile level of the 
distribution of limiting masses. These fits can be taken as esti- 
mates of the minimum masses. In Table |2] we list the values of 
these fits for different redshift bins and they will be considered 
in the following as miminum masses for different levels of com- 
pleteness. 

3.3. Star formation rates 

For the parent SDSS sample, star formation rates were taken 
from the MPA/JHU catalogue. These SFRs were derived us- 
ing the technique described in Brichmann et al. (2004). For the 



selected star-forming galaxies of the zCOSMOS sample, SFRs 
were calculated using the luminosity of the brightest available 
Balmer emission line. These values were corrected for aperture 
effects using the factors derived from HST photometry. In those 
cases in which this value was not available, we used Subaru pho- 
tometry. Extinction correction was carried out using the Balmer 
decrement for those objects with more than one Balmer hydro- 
gen recombination line with good S/N, and assuming the theo- 
retical ratios at standard conditions of temperature and density 
from Storey & Hummer (1995) and the CardelH et al. (1989) ex- 
tinction law. For objects with only one available Balmer line (a 
36% of the selected sample), we assumed a reddening coefficient 
from the E(B-V) parameter derived from the stellar synthesis fit- 
ting. The discrepancies found in those objects with both sources 
of information about the inner extincion do not vary the SFR 
distribution significantly. However, it must be stressed that ex- 
tinction of the gas and stars do not have to correlate, as is the 
case in this sample, so we preferred to use gas extinction when- 
ever possible. 

We derived a linear relation between Ha luminosity and the 
log(SFR) calculated by Brinchmann et al. (2004) for the DR7- 
SDSS sample in the Local Volume. This relation yields 

log(SFR) = 0.982 • log L(Ha) - 40.66 (4) 

and then we derived SFRs in the star-forming selected zCOS- 
MOS galaxies. In those redshift ranges (i.e. for z > 0.44) for 
which Ha is not observed in the VIMOS spectra, we used Hfi, 
Hy, or H<5 emission lines instead, based on the theoretical coef- 
ficients between Balmer hydrogen lines from Storey & Hummer 
(1995) for typical values of both electron temperature and den- 
sity. 

In Fig. [3] we compare the SFR distribution of the SDSS 
galaxies with the corresponding SFR distributions for different 
redshift bins in the zCOSMOS selected sample. The vertical blue 
solid line shows the median SFR for the whole zCOSMOS sam- 
ple (0.75 Mo/yr), while the red vertical solid line shows the me- 
dian SFR for each plotted bin. As can be observed, on one hand 
the median value for the whole zCOSMOS sample is higher than 
for the SDSS (with a median of 0.25 Mo/yr) and the median SFR 
in each bin increases with z from a value log(SFR) = -0.41 Mg/yr 
in the range < z < 0.2 up to log(SFR) = 1 .93 Mo/yr for z > 1 .0. 
This increase in SFR with cosmic time has already been estab- 
lished well by several authors (e.g. Daddi et al., 2007; Bouche et 
al., 2010), but may be also a selection effect when observing the 
brightest galaxies at larger z. 

3.4. Metallicity 

The metallicity of the studied star-forming galaxies can be es- 
timated using the oxygen abundances of their gas phase as a 
proxy. The most accurate method of deriving this chemical abun- 
dance is based on the previous determination of the electron tem- 
perature of the gas, via the intensity ratio of nebular-to-auroral 
emission lines (e.g. [Om] AA 5007, 4363) and the relative in- 
tensity of the strongest nebular emission lines to a hydrogen re- 
combination line (the so-called T^-method). However, as no in- 
formation is available in the VIMOS spectra of the zCOSMOS 
20k sample about the weak auroral lines, it is necessary to apply 
strong-line methods. These are based on direct calibration of the 
relative intensity of the strongest collisionally excited emission 
lines with grids of photoionisation models or samples of objects 
with an accurate determination of the oxygen abundance. 

In this work, our aim is to derive metallicities that agree 
with those derived using the Tg-method. Therefore, for the SDSS 
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Fig. 3. Distribution of the calculated SFR for the star-forming galaxies of the SDSS (in black) and for zCOSMOS star-forming 
galaxies in different redshift bins (in red). The vertical lines show the medians for the whole zCOSMOS sample (in blue) and for 
each redshift range (in red). 



le-K)5 



10000 !- 



100 r 



le-K): 



^ 10000 r 



E 




1000 r 



100 



7 8 9 

12-l-log(0/H) 



12-l-log(0/H) 



12-^log(0/H) 



Fig. 4. Distribution of the calculated 12-l-log(0/H) for the star-forming galaxies of the SDSS (black histogram) and zCOSMOS 
(red histograms) for different redshift bins. The vertical lines show the medians for the whole zCOSMOS sample (in blue) and 
for each redshift range (in red). 



sample and the zCOSMOS galaxies in the redshift range 0.01 < 
z < 0.44, we used the calibration of the N2 parameter proposed 
by PMC09, which is based on a sample of objects with a good 
determination of 12-i-log(0/H) via the Tp-method. The N2 pa- 
rameter, defined as the ratio of [Nii] A 6584 to Ha by Storchi- 
Bergmann et al. (1999), has been already used by Denicolo et al. 
(2002) as a metallicity proxy. The main advantages of this pa- 
rameter are i) its independence on reddening and flux calibration 
and ii) its monotonical relation with Z over a very wide range of 
metallicity. The PMC09 calibration, as discussed in Queyrel et 
al. (2012), is also consistent with other calibrations of the N2 
parameter based on photoionisation models (e.g. Denicolo et al., 
2002; Pettini & Pagel, 2004; Nagao et al., 2006). According 
to Perez-Montero & Diaz (2005), the main drawbacks of this 
parameter are its dependence on the ionisation parameter, the 
equivalent eff'ective temperature, and the N/O. Taking all these 
effects into account, the N2 parameter presents an overall uncer- 
tainty of about 0.3 dex across the entire metallicity range. This 
makes this parameter to be very uncertain for the determination 



of chemical abundances in single objects, but it offers many ad- 
vantages when it is used with statistical purposes. 

For higher redshifts, as the [Nii] emission-line is not acces- 
sible anymore in VIMOS spectra, we derived metallicities using 
the R23 parameter This was defined by Pagel et al. (1979) as 
the relative sum of [0 11] A Till and [O iii] AA 4959, 5007 to H/? 
intensities. Its relation with Z is bivaluated {i.e. it increases with 
increasing Z at low Z and it decreases with increasing Z at high 
Z) and it has a strong additional dependence on ionisation pa- 
rameter and equivalent effective temperature. To minimise this 
dependence we used the calibration proposed by Kobulnicky et 
al. (2003) based on the photoionisation models from McGaugh 
(1991). This calibration takes into account the dependence of 
R23 on the ionisation parameter with additional terms as a func- 
tion of the [O ii]/[0 in] ratio. The choice of the calibration branch 
(low- or high-Z) for each object was done according to the value 
providing less dispersion in the resulting mass-metallicity rela- 
tion. Using this criterion, we selected the upper-branch calibra- 
tion in all cases translating into oxygen abundances with a lower 
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limit of 12+log(0/H) = 8.0. However, it is expected that a sub- 
sample of the objects have lower metallicities, as those identified 
by Maier et al. (in prep.) using ISAAC-VLT near-IR observa- 
tions of zCOSMOS galaxies of the bright sample and measuring 
emission-lines which allow the choice of the appropriate branch 
of the R23 calibration. An estimate of the relative weight of these 
metal-poor objects in the z bins where double-valued strong-line 
methods were used can be obtained from the local SDSS sample. 
The proportion of objects with 12-i-log(0/H) < 8.0 from the N2 
method is only 0.2%, and if we take into account the minimum 
mass at 50% of completeness for z > 0.5, is only a 0.016%. 

Regarding the reddening dependence for the [O ii]/HyS ratio, 
we used the same procedure based on equivalent widths and the 
factor or as a function of D(4000), as described in section [TT] 

Finally, for galaxies with a redshift z > 0.9, we used the 
method described by Perez-Montero et al. (2007) based on the 
relation between [Neiii] A 3869 and [Om] A 5007 (the so-called 
02Ne3 parameter) and the theoretical relation between H/3 and 
other Balmer lines at bluer wavelengths with enough S/N (i.e. 
Hy and H6). For these objects, we also considered reddening 
correction from the Balmer decrement when more than a hydro- 
gen emission-line is available and, otherwise, we took the stellar 
extinction derived from spectral synthesis fitting as described in 
section im 

As already described by Kewley & Ellison (2008), impor- 
tant differences can arise between the metallicities derived using 
different strong-line methods and/or different calibrations of the 
strong emission-line ratios. In our case, as the PMC09 calibra- 
tion of the N2 parameter is based on the compilation of objects 
with a "direct" determination of the oxygen abundance, while 
the calibration of R23 is based on sequences of photoionisation 
models, this difference can be very large. To convert the metallic- 
ities derived from R23 and 02Afe3 parameters to those estimated 
in the local universe from the N2 parameter, we used the fol- 
lowing linear relations which is based on models by Chariot & 
Longhetti (2001) as described in Lamareille (2007). 

log(0/H)N2 = 1 .228 X log(0/H)R23 - 1 .92 1 (5) 

for the lower branch calibration of R23 and 

log(0/H)N2 = 1 .203 X log(0/H)R23 - 1 .932 (6) 

for the calibration of the upper branch. 

In Fig.|4]we show the oxygen abundance distribution for the 
SDSS selected star-forming galaxies compared with the same 
distributions for different z bins in the zCOSMOS sample in or- 
der to investigate their completeness in metallicity. In the low- 
est redshift bins (z < 0.4), the high-metallicity values measured 
in SDSS galaxies are not reached in zCOSMOS basically be- 
cause almost no massive galaxies are detected, as the sampled 
volume probed by zCOSMOS at these low redshifts is much 
smaller than for SDSS. For low metallicities, no values lower 
than 12-t-log(0/H) < 7.7 are found in zCOSMOS in the same 
low z bins because the minimum [N 11] flux which can be mea- 
sured in this survey is brighter than in SDSS. Nonetheless, the 
statistical weight of this low-metallicity queue in the SDSS sam- 
ple is negligible (0.02%). We can thus consider that our selection 
of zCOSMOS galaxies is representative, in terms of metallicity, 
of the low-redshift (z < 0.4) star-forming galaxies. 

At larger redshift (z > 0.4), there is also a selection effect re- 
lated to the minimum emission-line flux detected in zCOSMOS, 
but this time it affects the high-metallicity regime. Indeed, in 
this regime, all oxygen abundances were derived from the upper- 
branch calibration of R23 (or 02Ne3 for z > 0.9), where weaker 



relative emission lines imply higher Z. Therefore, the minimum 
detectable [0 11] and [O iii] emission lines imposes an upper limit 
on the metallicity of zCOSMOS galaxies. However, as in the 
case of low-metallicity galaxies at low redshift, this has a neg- 
ligible impact on the final distributions as SDSS galaxies with 
12+log(0/H) > 9 only represents 1.14 % of the total sample. 
Regarding the cutoff seen for lower metallicity at this redshift 
regime, it corresponds to the choice of the upper branch cali- 
bration in all cases (12-i-log(0/H) > 8.0). Although a number of 
galaxies is expected to have low metallicities (as those observed 
with ISAAC in Maier et al. (in prep.)), these galaxies have a 
negligible statistical weight when we compare with the SDSS 
sample (0.2%). In any case, this number cannot be taken as an 
estimate of the relative number of very-low metalUcity galaxies 
in all the zCOSMOS sample. 

3.5. Nitrogen-to-oxygen abundance ratio 

The N/O was derived for both the SDSS sample and the zCOS- 
MOS 20k sample (up to z ~ 0.42 only), using the empirical 
calibration of the N2S2 parameter. This parameter is defined as 
the ratio of the emission-line fluxes of [N 11] A 6584 and [S 11] AA 
6717, 6731, which according to PMC09 has a linear relation with 
log(N/0). Although the ratio between [Nii] A 6584 and [On] A 
3727, the so-called N202 parameter, is also a good indicator of 
N/O (Perez-Montero & Diaz, 2005; PMC09), we did not use it 
for sake of consistency with the zCOSMOS data, for which no 
[0 11] information exists for z < 0.5. 

In Fig. [5] we show the distribution of log(N/0) for both the 
SDSS and zCOSMOS samples. Two redshift bins are considered 
for the zCOSMOS sample: z < 0.2 and 0.2 < z < 0.42. Contraiy 
to metallicity, the N/O is well sampled for low values in zCOS- 
MOS using the N2S2 parameter. However, this is not the case 
for high values of N/O. As in the case of metallicity, this is be- 
cause only few massive galaxies are sampled at these redshifts 
within zCOSMOS, and thus very high values of N/O are missing 
in this sample. This effect is also seen in the zCOSMOS sample 
because the median log(N/0) for the lower redshift bin (-1.38) 
is lower than for z larger than 0.2 (-1.29), and both of them are 
much lower than the median log(N/0) for the entire SDSS sam- 
ple (-0.83) since more massive galaxies are sampled on it. 

4. Results and discussion 

4.1. Cosmic evolution of fhe mass-metallicity relation 

To consistently study the cosmic evolution of the MZR we de- 
rived its shape in the selected star-forming galaxies of the SDSS, 
which are representative of the local universe. As described in 
the sections above we used the compiled stellar masses from the 
MPA/JHU catalogue and calculated oxygen abundances from 
the PMC09 calibration of the N2 parameter The derived MZR 
for the SDSS is shown in Fig. |6] We computed a quadratical fit 
to the median values of O/H per stellar mass bins of 0.2 dex. The 
resulting fit has the following expression 

3/ = 1.1598 + 1.2971 ■x-0.0545 -jc^ (7) 

where y is 12-i-log(0/H) and x is log(M») in units of M©. This 
fit is vaHd in the range of steflar mass 8.0 < log(M/MQ) < 1 1.8. 
The dispersion was calculated as the standard deviation of the 
residuals to the resulting fit, giving a result of 0.090 dex. This 
dispersion is slightly larger for log(M») < 10 (0. 105 dex) than for 
higher stellar masses (0.081 dex). We checked that our quadrat- 
ical fit is consistent with the results given by Kewley & Ellison 
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Fig. 5. Distribution of the calculated log(N/0) for the star-forming galaxies of the SDSS (black histogram) and zCOSMOS (red 
histograms) for different redshift bins. The vertical lines show the medians for the whole zCOSMOS sample (in blue) and for 
each redshift range (in red). 
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Fig. 6. Relation between stellar mass (in units of M©) and 
12+log(0/H) for the selected star-forming galaxies of the SDSS- 
DR7. The red solid line shows the quadratical fit to the medians 
in bins of 0.2 dex of stellar mass. 



(2008) for a similar subsample of the SDSS. By converting our 
oxygen abundances to those studied in that work using different 
strong-line methods or calibrations of N2, we measured that our 
median Z is 0.10 dex higher on average for all stellar masses. 
However, we also found that by increasing the S/N cutoff of our 
sample to 10, the agreement with the fittings of the MZR pro- 
vided by Kewley & Ellison (2008) is perfect. Therefore, we con- 
clude that using a high S/N cutoff introduces a selection effect in 
the subsample, which mainly affects to the metal-rich galaxies 
when it is applied to all the involved emission lines. 

In Fig. |7] we show the MZR for six different redshift bins 
between 0.01 and 1.32 and compared to the local MZR derived 
above. Each panel also shows the ±cr curves of the medians of Z 
in each mass bin and the minimum masses for different levels of 
completeness as derived in section 3.2. In the three panels corre- 
spondinbg to z between 0.4 and 1 .0 we show the small subsam- 



ple of zCOSMOS galaxies observed in the near-infrared (Maier 
et al., in prep.) and for which oxygen abundances from the lower 
calibration of the R23 parameter have been derived. 



In Table [3] we summarise for each redshift bin the vertical 
offset between the local MZR and the average of medians for 
different mass intervals in each z bin. We also give the corre- 
sponding dispersions, calculated as the standard deviation of all 
the medians. In all cases, for the sake of consistency, these offsets 
were calculated taking from the SDSS those subsamples of ob- 
jects in the same metallicity and stellar mass range measureable 
in each z bin of zCOSMOS. The inspection of the values listed 
in that table leads to the following, (i) The derived median oxy- 
gen abundance for the same bins of stellar mass increases with 
cosmological age so that it was around half the present value at 
z a! 1 ( » 8 Gyr ago), (ii) The increase in the oxygen abundance 
with cosmic time is not uniform but, on the contrary, it looks to 
be significantly larger above z ~ 0.5. (iii) Both the dispersion 
associated with the derived medians, and the oxygen abundance 
uncertainties intrinsic to the strong-line methods for each z bin, 
are larger than the offset for z < 1.0. 



Some authors {e.g. Kobulnicky et al. 1999) have warned 
that taking too low a threshold of S/N for the Balmer emission 
lines introduces significant uncertainties in the determination of 
metallicity using strong-line methods, such as R23. On the con- 
trary, as explained in previous sections, taking too high a thresh- 
old for certain lines can introduce a bias in the selected sample. 
To quantify how changing the S/N threshold for Balmer lines 
affect the median Z in the different redshift bins, we calculated 
them taking only Balmer emission-lines with S/N > 5, instead 
of 2. Our results indicate that the number of selected galaxies is 
considerably decreased for the high z bins (more than 60% for z 
> 1 and more than 35% for z > 0.4). In Table[3]we also list the Z 
offsets in relation to the SDSS MZR. All of them give higher Z 
than in relation to the samples with S/N threshold of two for all 
the involved emission lines, and this difference increases with z, 
as an increasing number of low-Z galaxies is excluded in each 
bin. However, the trend to finding lower median Z with higher 
z is maintained, although serious caveats must be taken for the 
intensity of this evolution when selection effects due to different 
S/N threshold criteria are imposed. 



9 



E. Perez-Montero et al.: Cosmic evolution of galaxy metal content over the last 10 Gyrs 

Table 3. Average offsets between the medians measured in the SDSS reference sample and for the star-forming galaxies of 
zCOSMOS 20k in each redshift bin for the mass-metallicity relation (MZR), the star formation-corrected mass-metallicity rela- 
tion (SMZ), and the relation between Stellas mass and N/O (MNOR). The dispersions of the fits are also given. 



z bin 


A MZR 


a MZR 


AMZl^ 


ASMZR 


crSMZ 


A MNOR 


a MNOR 


SDSS 


0.00 


0.09 




0.00 


0.08 


0.00 


0.14 


0.0 - 0.2 


-1-0.01 


0.16 


+0.01 


+0.01 


0.15 


0.00 


0.23 


0.2 - 0.4 


-0.04 


0.20 


-0.02 


-0.01 


0.19 


-O.lf] 


0.31" 


0.4 - 0.6 


-0.14 


0.19 


-0.11 


-0.12 


0.19 






0.6 - 0.8 


-0.19 


0.17 


-0.12 


-0.13 


0.17 






0.8- 1.0 


-0.17 


0.19 


-0.08 


-0.10 


0.19 






1.0-1.32 


-0.36 


0.15 


-0.27 


-0.30 


0.15 







" with B aimer lines with S/N > 5 
* up to z a 0.42 




log(M,/M ) log(M,/M ) log(M./M ) 

Fig. 7. Cosmic evolution of the MZR for six different redshift bins increasing from left to right and from up to down. The red 
solid line shows the fit to the MZR in the SDSS data. Grey areas show the +cr intervals to the Z medians fits in different mass 
bins for each redshift range. The vertical lines show the minimum masses for 25%, 50% and 75% levels of completeness. Red 
crosses represent low-Z objects in the sample found using ISAAC near-lR observations (Maier et al. in prep.) 



4.2. The SFR-corrected mass-metallicity relation and its 
evolution 

Stellar mass and metallicity are two of the main fundamental 
parameters of galaxies. Their relationship for complete samples 
of galaxies, as shown in the above section, gives us valuable in- 
formation about the main mechanisms governing the evolution 
of galaxies. However, recent works {e.g. Mannucci et al., 2010; 
Lara-Lopez et al. 2010) have demonstrated the important role of 
SFR, which is also tightly related to stellar mass and metallicity 
and thus affects the shape and the dispersion of the MZR. For 



instance, galaxies with higher SFR have systematically lower 
metallicities. This effect could be explained by invoking the in- 
teraction of galaxies with inflows of unenriched gas that trig- 
gers the episodes of star formation and reduces the relative abun- 
dance of oxygen in the gas phase. For this reason, Mannucci et 
al. (2010) defines a fundamental metallicity relation (FMR) that 
is basically a surface in the three-dimensional space formed by 
stellar mass, Z, and SFR in which neither selection effects due 
to SFR appear nor evolution with z are detected up to z ~ 2.5. 
This lack of evolution in the defined 3D surface (although other 
authors, such as Lara-Lopez et al. (2010) defined it as a plane) is 
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Fig. 8. To the left MZR and quadratical fits for different bins of log SFR. To the right, same plot for the relation between log 
(N/O) and stellai^ mass. 
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Fig. 9. Relation between stellar mass in units of M© and 
12+log(0/H) with the correction depending on log(SFR) for the 
star-forming selected galaxies of the SDSS-DR7. The red solid 
line shows the quadratical fit to the medians in bins of 0.2 dex of 
stellar mass. 



also found for a subsample of the zCOSMOS-bright galaxies in 
the redshift range 0.2 < z < 0.8 (Cresci et al. 2012). However, in 
this paper, since we try to isolate the metallicity evolution with 
cosmic age independently of the stellar mass, we redefine the 
question, isolating the dependence of Z with SFR to remove the 
selection effects at high z- In the left hand panel of Fig. [8] we 
show the MZR with different quadratical fits for different bins 
of SFR for the aforementioned selected star-forming galaxies of 
the SDSS. As mentioned above, there is a trend toward obtaining 



lower 12H-log(0/H) for higher SFRs. We then applied the SFR 
correction only to the metallicity, in order to remove possible se- 
lection effects at higher z- We thus used the SFR-corrected MZR 
(hereafter S MZ relation). 

By inspecting Fig.[8]again, we clearly see that the evolution 
of the vertical offsets between the quadratical fits to subsamples 
with different SFRs is not linear Indeed, we can fit a quadratical 
SFR-correction to the Z axis in order to reduce the dispersion 
and remove the dependence on SFR. The additional term to add 
to the metallicity axis is 



/(log SFR) = 0.0505 ■ (log SFR) - 0.0144 ■ (log SFRf. 



(8) 



In Fig.|9]we show the corrected MZR, in which the corrected 
metallicity is represented as a function of the stellar mass. The 
resulting quadratical fit to the Z medians in bins of 0.2 dex in 
stellar mass in the range 8.0 < log(M/MQ < 11.8 is 



y = 0.4655 h- 1.41 18 ■ x - 0.0589 ■ 



(9) 



where y — 12H-log(0/H)-i-/(SFR) and x is the stellar mass in units 
of solar masses. The dispersion of this new relation is slightly 
lower than for the uncorrected MZR (0.083 dex) and, as in the 
uncorrected case, the dispersion is lower for masses higher than 
log(M,) = 10.0 (0.075 dex) than for lower masses (0.097dex). 
An explanation of the fact that the dispersion is only reduced 
by 0.01 dex at all stellar masses can be found the majority of 
SDSS galaxies not having extreme SFR values. The dispersion 
in the SFR distribution of the SDSS sample is only 0.37 dex. 
A much more thorough calculation of this correction should be 
done using well-known galaxies with extreme values of the SFR 
to put more precise constraints on the SFR-correction term of 
the 5 MZ. 

To study to what extent selection effects of the SFR can af- 
fect the conclusions reached from the study of the MZR, we ap- 
plied the SFR-correction to the selected zCOSMOS galaxies at 
different z bins. The SMZ is plotted in Fig. [TO] for the same six 
different redshift bins between 0.01 and 1.32 with respect to the 
relation derived for the SDSS sample. The panels also show the 
±lcr bands around the fitting to the medians for Z-i-/(SFR) in 
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Fig. 10. Cosmic evolution of the SFR-corrected MZR for six different redshift bins increasing from left to right and from up to 
down. The red solid line shows the fitting to the S MZ in the SDSS data. Upper and lower solid lines show the ±cr intervals 
to the fitting to the y medians in different mass bins for each redshift range. The vertical lines show the minimum mass limits 
for 25%, 50% and 75% levels of completeness. Red crosses represent low-Z objects in the sample found using ISAAC near-IR 
observations (Maier et al. in prep.) 



each mass bin and the minimum mass for different levels of com- 
pleteness. In Table[3] we summarise for each redshift bin the ver- 
tical offset between the reference SDSS S MZ and the average of 
median offsets for each mass bin and the dispersion, calculated 
as the standard deviation of the residuals to the fitting. 

The inspection of the values listed in that table leads to the 
following conclusions, (i) There is still a trend toward lower me- 
dian values of the SFR-corrected Z for higher z, although the 
measured offsets are lower This implies that, although the MZR 
is affected by SFR selection effects for samples taken at high 
z, these are not sufficient to explain the observed evolution, as 
stated in other studies (e.g. Mannucci et al. 2010; Cresci et al. 
2012). (ii) The dispersion of the SFR-corrected Z is not notice- 
ably lower in each z bin than for the MZR, possibly because that 
the SFR distributions for each z bin do not have very large dis- 
persions, so the correction affects the median value but not the 
dispersion of each distribution. 

4.3. The relation between stellar mass and 
nitrogen-to-oxygen ratio and its evolution 

The N/0 complements the study of the chemical composition 
of the nebular phase of star-forming galaxies, and it is a robust 
indicator of the chemical evolution of galaxies. Besides, as the 
production rate of secondary N is faster than for O (Henry et al. 
2000; Thuan et al. 2010), the amplitude of the evolution of its re- 



lation with stellar mass gives more reliability to the conclusions 
reached about its evolution with cosmic age. 

Although the relation between stellar mass and N/O (here- 
after MNOR) has been already studied by PMC09, we re- 
calculated this relation using the same selected star-forming 
galaxies of the SDSS in order to be compared with the zCOS- 
MOS galaxies at higher z. In Fig. 11 we show this relation be- 
tween stellar mass and log(N/0) for the selected star-forming 
galaxies of the SDSS-DR7. As explained above, the N/O was 
calculated using the PMC09 calibration of the N2S2 parameter. 
As in the case of the MZR, there is a trend toward finding higher 
N/O for more massive galaxies because all galaxies in this sam- 
ple akeady lie in the regime of production of secondary N and 
so its N/O increases with metallicity. We performed a quadrati- 
cal fit to the medians of N/O in bins of stellar mass of 0.2 dex in 
the stellar mass range 8.6 < log(M/M0) < 11.6 which gives the 
following expresion: 



y = -9.4457 + 1.3577 ■ x - 0.0501 ■ 



(10) 



where y is log(N/0) and x is log(M») in units of solar masses. 
The dispersion, calculated as the standard deviation of the resid- 
uals, is equal to 0.144 dex. Contrary to the MZR relation, this 
dispersion is slightly lower for masses greater than 10'° solar 
masses (0.136 dex) than for lower masses (0.156 dex). The more 
restricted stellar mass range where the fit is possible and the 
higher dispersion for lower stellar masses are possibly owing to 
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Fig. 11. Relation between stellar mass in units of Mg and 
log(N/0) for the star-forming selected galaxies of the SDSS- 
DR7. The red solid line shows the quadratical fit to the medians 
in bins of 0.2 dex of stellar mass. 



the presence of objects with additional primary N production at 
the low stellar mass regime. 

Besides, unlike the MZR too, there is no apparent selection 
effect due to different values of the SFR, as can be seen in right 
hand panel of Fig. [8] This lack of SFR-dependence of the rela- 
tion between stellar mass and N/O could be a confirmation that 
the relation between Z and SFR is mainly due to the inflows of 
metal-poor gas, which would trigger the star formation, decreas- 
ing the relative abundance of the metals but which, in contrast, 
do not affect the ratio between secondary and primary elements 
(Edmunds 1990). 

Therefore, the study of the evolution of the MNOR with cos- 
mic time represents a robust tool for studying the evolution of 
the chemical content of galaxies as a function of the age of the 
Universe, as N/O does not depend on selection effects related to 
SFR. In Fig. [12] we show the evolution of this relation for the 
two redshift bins in zCOSMOS for which we were able to de- 
rive the N/O using the N2S2 parameter The offsets between the 
relation derived in the SDSS and for the two analysed z bins and 
their dispersions are listed in Table [3] The analysis of the ob- 
tained values leads to the following conclusions, (i) For a given 
stellar mass, there is a significant evolution in the median N/O 
abundance ratio of galaxies from z ~ 0.4 up to now. (ii) This 
evolution is a consequence of the evolution of the metallicity of 
galaxies with cosmic age. As can be seen in Fig. 13 the aver- 



age relation between 12-i-log(0/H) and log(N/0) in the regime 
of production of secondary nitrogen is kept in the two redshift 
bins, (iii) This evolution cannot be due to any selection effect 
of the sampled galaxies at high redshift that have, on average, 
larger SFR, because N/O is not affected by inflows of metal- 
poor gas in galaxies, (iv) In the secondary nitrogen production 
regime, N/O increases linearly with oxygen abundance, as can 
be seen in Fig. 13 Therefore the N/O evolution derived from z ~ 
0.3 can be transformed into metallicity evolution using the slope 
of this linear relation. This slope is around 1 .4 for both the two 



redshift bins analysed, so this gives a variation in 12H-log(0/H) 
of 0. 1 dex on average, which is slightly higher than the values 
found using both MZR and SMZ. However, it must be taken the 
large scatter into account in the N/O vs. O/H when computing 
the O/H evolution in this way 



5. Summary and conclusions 

We studied the cosmic evolution of the MZR from z ~ 1.3 up 
to now using complete samples of star-forming galaxies: in the 
Local Universe we used the SDSS-DR7 and for higher z we used 
the 20k zCOSMOS bright sample. The selection of star-forming 
objects and derivation of physical properties, including stellar 
masses, star formation rates, oxygen abundances, and N/O, were 
performed using consistent methods. This is especially relevant 
in the case of the VIMOS spectra taken for the zCOSMOS sam- 
ple since the set of available emission lines varies in each redshift 
regime. Two important aspects were considered when deriving 
chemical abundances: i) all strong-line methods used to derive 
12H-log(0/H) and log(N/0) were consistent with the T^-method 
as used in PMC09, and ii) the S/N threshold to select objects was 
put at a lower value to avoid possible selection effects in the final 
distributions of the derived properties. 

The analysis of the evolution of the MZR in different redshift 
bins confirms an increase in the metallicity from z-l .32, but the 
dispersion is too high to assess any variation in the shape of the 
MZR. 

Consideringthe dependence of metallicity on SFR found by 
several authors, we derived an SFR-corrected MZR. This cor- 
rection slighly reduces the dispersion of the MZR, but not no- 
ticeably since the dispersion in the SFR distributions is not very 
high. We still see the evolution of the metallicity for galaxies 
of the same stellar mass bin at different epochs, although this is 
slightly lower than for the MZR. 

Finally, we studied the evolution of the MNOR taking the 
log(N/0) instead of (O/H) in the MZR. Since N/O was derived 
from the N2S2 parameter as calibrated by PMC09, we carried 
out this study for the zCOSMOS 20k sample only up to z ^ 
0.42. The main advantage to considering the MNOR with re- 
spect to the MZR is the independence of this relation on SFR, 
which confirms that the dependence of Z on SFR is probably 
due to inflows of metal-poor gas in some starburst events. In the 
studied redshift regime we still see the evolution of the chemical 
content of galaxies with the same stellar mass from younger ages 
of the Universe, confirming that the aforementioned evolution of 
the MZR is not just a selection effect of the observed objects at 
high z. 
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